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CHARTS AND APPROXIMATE FORMULAS FOR THE ESTIMATION OF AEROELASTIC EFFECTS 
ON THE LATERAL CONTROL OF SWEPT AND UNSWEPT WINGS' 

By Kenneth A. Foss and Franklin W. Diederich 


SUMMARY 

Charts and approximate Jormulas are presented for the estima- 
tion of static aeroelastic effects on the spanwise lift distribution, 
roliiiKj-moment coefficient, and rate of roll due to the deflection 
of ailerons on swept and unswept wings at subsonic and super- 
sonic speeds. Some design considerations brought out by the 
results of this report are discussed. 

This report treats the lateral-control case in a manner similar 
to that employed in NACA Report 1140 for the syinmetric- 
f light case and is intended to be used in conjunction with NACA 
Report 1140 and the charts and formulas presented therein. 

INTRODUCTION 

'riio laloral control and maneuvorabilit}' of a wing are 
important design considerations. "Hiese characteristics may 
he affected to a significant extent aeroelastic action, 
parliciilai’ly at liigli dynamic pressures and in the case of 
lliin wings, swept wings, and wings designed for low wing 
loadings, because the opei*ation of ailerons and spoilers 
usually creates aerodAuiamic forces which deform the wing. 

As a result of these deformations, the angles of attack 
along tile span often change in such a maimer as to pi'oduce 
lifts which oppose the rolling moment of the aileron or 
spoiler; furthermore, these lifts cause additional deforma- 
tions which may again reduce the rolling moment, and so on, 
until equilibrium is reached. Wing flexibility may thus 
cause a serious loss in the control power; in fact, if the 
dynamic jiressure of the airstream is sufficiently high, the 
aileron rolling moment may be completely nullified. The 
speed and dynamic pi-essure at this condition are often re- 
ferred to as the aileron reversal speed and reversal dynamic 
pressure, because at higher dynamic pressures the conti-ois 
would have to be reversed in order to roll the airplane. 
Wlien wing flexibility causes a loss in lateral control, there is 
also usually a loss in the rolling maneuverability, which may 
be expressed as the wing-tip helix angle due to rolling and is 
affected by changes in both the control power and the damp- 
ing in roll. These aeroelastic efi'ects on the lateral control 
and maneuverability have to be taken into account in the 
design of a wing. 

Several methods are available for calculating these 
effects (ref. 1, for instance), but, since these effects depend 


on the structural characteristics of the wing, wliich are not 
accurately Imovm in advance of its design, the relatively- 
large amount of time required for even the most efficient 
of these methods militates against theii* use in coiuiection 
with preliminary design calculations. A need exists, there- 
fore, for means of estimating quickly some of the more 
important aeroelastic effects on lateral control with an 
accuracy' that is sufficient for preliminary design pinposes. 

The related problem of estimating static aeroelastic 
eft'ects on the magnitude and spanwise distribution of the 
lift in symmetric flight has been treated by' the charts and 
approximate formulas presented in reference 2. Tlie pres- 
ent report consists of an extension of the anah'sis of refer- 
ence 2 to the lateral-control case. Inasmuch as the static 
aeroelastic equations are linear, the results presented in 
the two reports may* be superimposed. Included in tlie 
present report are approximate formulas for the estimation 
of the static aeroelastic effects on the spanwise lift distri- 
bution, rolling moment, and rate of roll due to aileron deflec- 
tions on swept and unswept wings at subsonic and supersonic 
speeds. Also presented are summary' charts which 
indicate whether a given design is likely’ to be afl*ected by 
losses in lateral control. By' means of these charts and 
approximate formulas as well as those of reference 2, the 
conventional procedure of designing a wmg on the basis of 
certain strength criteria, checking it for aeroelastic phe- 
noiTiena, and tlien reinforcing it, when necessary', to meet 
the stiffness requirements imposed by' these phenomena can 
often be simplified greatly', inasmuch as the efiect of some of 
these phenomena can be estimated m advance of design. 

In order to keep the length of the report to a minimum 
and to avoid a repetition of much of the material presented 
in reference 2, the present report has been ml t ten in such a 
manner as to facilitate its joint use with reference 2 ratlier 
than to make it entirely’ self-contained. The use of the charts 
and approximate formulas presented herein is described and 
the limitations of the charts and the light they' shed on some 
design problems are discussed. A numerical example is 
included to illustrate the use of the approximate formulas 
of this report. A brief description of the calculations on 
which the charts and approximate formulas are based is 
contained in the appendix to supplement the more detailed 
derivations in references 1 and 2. 


• Previously released as NACA TN 27*17, “Charts and .Approximate Formulas for the Estimation of Aeroelastic Effects on the Lateral Control of Swept and Unswept Wings*' by 
Kunneih A. Foss and Franklin W. Diederich, 1952. 
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SYMBOLS 


aspect ratio, b^jS 

swept-span aspect ratio, ^1/cos-A 
location of section aerodynamic center 
measured from leading edge, fraction 
of cliord 
wing span, in, 

chord (measured perpendicular to elastic 
axis), in. 

aileron chord, in. 
section lift-curve slope per radian 
effective wing lift -curve slope per radian 
rolling-moment coefficient 
damping-in-roll derivative 
rolling-moment coefficient due to aileron 
deflection 

location of chordwise center of pi*essurc 
of lift produced Idv aileron deflection 
measured from leading edge, fraction 
of chord 

dimensionless sweep parameter. 


{GJ)r 


tan- A 


{EI)r 

Young’s modulus of elasticity, Ib/sq in. 
location of elastic axis measured from 
leading edge, fraction of chord 
dimensionless moment arm of section lift 
about elastic axis, e—a 
dimensionless moment arm of lift due to 
aileron deflection about elastic axis, 


cpi—e 

root -St iff ness function given in equation 
(B25) of reference 2 
allowable bending stress, Ib/sq in. 
dimensionless function of distance along 
span used in approximate formulas for 
angle of attack due to aeroelastic 
action of ailerons 
modulus of rigidity, Ib/sq in. 
wing thickness, in. 

section bending moment of inertia, in.** 
mass moment of inertia of entire aii-plane 
about its longitudinal axis, in.^ 
mass moment of inertia of both wings 
about longitudinal axis of airplane, in.*^ 
section twisting moment of inertia, in.^ 
dimensionless parameters used in approx- 
imate formulas for dimensionless dy- 
namic pressures at aileron reversal 
given in table I 

dimensionless sweep parameter, 

A; the dimensionless 

C,(V {EI)r 

parameter kje is identical to k except 
that (\ is replaced by 
lift per unit distance along s])an, Ib/in. 
rolling moment on botli wings, in-lb 
bending moment al)Out an axis perpen- 
dicular to elastic axis, in-lb 
free-stream Mach number 
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Subscripts: 


T) 


0 

H 


r 


t 

0 


design load factor 
rolling acceleration, radians/sec- 
wing-tip helix angle due to roll, radians 
(hmamic pressure, Ib/sq ft 
dimensionless d^mamic pressure, 

— i dimensionless 

144 {GJ)r 

dynamic pressure is identical to q* 
except that Ci is replaced b}’* co 
dimensionless dynamic pressure, 
g CL^CrSt^ sin A 

l44 {EI)r 
total wing area, sq in. 
distance along elastic axis measured from 
wing root, in. 

dimensionless distaiice along elastic axis, 
s/St 

accumulated torque about clastic axis, 
in-lb 

distributed torque due to inertia loading, 
in-lb/in. 
airspeed, ft/sec 

design gross weigliL of airplane, lb 
weight of primarv structure of both wings, 
lb 


lateral coordinate, in. 
angle of attack in planes parallel to plane 
of symmetry, radians 
angle of attack equivalent to unit aileron 
deflection, radians 

angle of local dihedral, radians; or span- 
wise slope of normal displacement of 
elastic axis 

aileron deflection measured in planes 
parallel to ah*stream, radians 
moment-arm ratio, e^lei 
structural-eftectivcncss factors defined in 
equation (15) of reference 2 
structural-efl'ectiveness factors defined in 
table 1 of reference 2 
angle of sweepback at clastic axis 
wing taper ratio, CtjCr 
density, slugs/cu ft 

angle of structural twist in planes per- 
pendicular to clastic axis, radians 

(Elh 


tip stiffness ratio. 




* con»tant f|re«« 

unit step function of distance along span 


at divergence 

at inboard end of aileron; inertia, in 
equation (36b) 

at outboard end of aileron; except in AIo 
at aileron reversal 
at wing root 

structural (due to structural deformation) 
at wing tip 

rigid wing (for €q*=q = 0) 
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USE OF THE CHARTS AND APPROXIMATE FORMULAS 

SUMMAUY OF METHOD AND SCOPE OF THE CALCULATIONS ON WHICH 
THE CHARTS AND APPROXIMATE FORMULAS ARE BASED 

A brief description of tlie method and scope of the cal- 
culations is given here lo indicate the limitations of the 
charts and approximate formulas. A detailed description 
of the method is given in the appendix. 

Most of these calculations were performed by an exten- 
sion, based on reference 1, of the matrix method described 
in appendix A of reference 2. This method consists in 
solving the dilferential equations descriptive of an elas- 
tically deforino<I wing under aerod\mamic loadings by 
numerical methods employing matrix technic|ues. Treated 
by this method were wings with three taper ratios (1, 0.5, 
0.2y one aileron configuration (50 percent semispan, out- 
board), two types of stifi’ness distributions, several values 
of the sweep parameters k and d which include sweptfor- 
ward, unswept, and sweptback wings, and several values of 
th(‘ section moment-arm ratio € and the dynamic-pressure 
ratio ([Iqn- Calculated for each case were tlie dynamic 
pressure at aileron reversal and the changes in the spanwise 
lift distribution and rolling moment due to aileron deflection. 
For the constant-choi-d — constant-stiffness wings, calcula- 
tions were also perfornuMl by an extension of the anal}Tic 
method described in appendix A of reference 2, which con- 
sists in solving the differential equations exactly for these 
relatively simple cases; these calculations were made for 
two aileron spans and several values of the parameters 
k, (ly €, and fn all cases the ratio of the aileron chord 

to the wing chord was constant along the span of the aileron. 

Some approximations liave been made in the calculations 
conc(‘rning the aerodpiamic induction effects, the root 
rotations, and the stiffness distribution, primarily in order 
to hold the number of variables considered in the analysis 
lo a minimum and to make the results more generally 
applicable. 

A(*rod 3 mamic induction effects at subsonic speeds are 
taken into account bv an overall reduction of the strip- 
theorv loading and, in the matrix calculations, b^' rounding 
off the strip-theory loading at the tip (see refs. 1 and 2); 
for supersonic speeds sti’ip theoiy is used with a small reduc- 
tion at the tip in the matrix calculations. This approxima- 
tion has made it unnecessary to consider e.xplicitl}' the 
effect. s of aspect ratio, sA\eep, and Madi number on the 
rigid-wing lift distribution; the effects of these parameters 
on (he total lift and on the aeroelastic increment to the lift 
distribution liave been taken into account. 

The rigid-bodj" rotations imparted to a SAvept wing bj' its 
triangular root portion vary among difl’erent designs in a 
largelA" unpredictable manner. They have therefore been 
taken into account only by the use of an effective root, the 
selection of which in any given case is discussed brieflA' in a 
subsequent section. 

'File spanwise distributions of the bending and torsional 
stiffnesses depend on the detailed design of the wing and 
cannot be generalized easily. The stiffness distributions 
used in the calculations of aeroelastic effects were obtained 
from tlie constant-stress concept outlined in appendix B of 
ref(‘rence 2, which constitutes an effort to relate the stiff- 


ness of a wing to its strength on the basis of the following 
assumptions: 

(1) The combined bending and torsional stresses are 
constant along tlie span. 

(2) The bending and torsional stresses are combined in 
such a manner that the sum of the ratio of the actual to the 
allowable bending stress and the ratio of the actual to the 
allowable torsion stress is equal to unity when the margin 
of safety is zero. 

(3) The structure is of the thin-skin, stringer-reinforced 
shell type and its main features do not vary along the span; 
for instance, the number of spars and their chordwise loca- 
tions are constant along the span. 

(4) At the design condition the spanwise distribution of 
the applied loading is proportional to the chord. 

Also used in the calculations Avere stiffness distributions 
which vary as the fourth power of the chord, as do those of 
solid wings and wings with geometrically similar cross 
sections; as pointed out in a subsequent section, the results 
of these calculations can be used to estimate aeroelastic 
phenomena of some wings which have large cutouts or 
which for some other reason do not liave stiffness distribu- 
tions represented fairly closely by those of tlie constant- 
stress typo. 

All calculations are based on the assumptions tliat twist- 
ing is resisted primarih' by the torsion cells of the wing 
structure and that the Aving deformations can be estimated 
by means of the elementary theories of bending and torsion 
about an elastic axis. 

SELECTION OF PARAMETERS 

Geometric parameters. — The geometric parameters used 
in the aual 3 ^sis are defined in figure 1. The location of the 
Center of pressure due to aileron deflection-.. 
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effective root indicated in tliis figure is discussed in refer- 
ence 2, In the present report the angle of aileron deflection 
8 is defined as being measured in planes parallel to the air- 
stream, This angle is equal to the product of the angle of 
rotation of the aileron about the hinge axis and tlie cosme of 
the sweep angle of the hinge axis. 

Although most of tlie charts and approximate formulas 
are based on a half-span outboard aileron (5i*=0.5, 5o*=l), 
the results of the analysis of the uniform wings with full- 
span ailerons (5f*=0, ^?o*=l) and “tip ailerons” 

5o*=l) indicate that, except for the angle-of-attack distri- 
butions, the results based on a half-span aileron may be 
expected to be valid for outboard ailerons having spans 
which differ considerably from one-half. 

Aerodynamic parameters. — The aerodynamic parameters 
which enter the analysis are the effective wing lift-curve 
slope, the location of the wing aerodynamic center, the loca- 
tion of the chord wise center of pressure due to aileron 
deflection, and the angle of attack equivalent to unit aileron 
deflection. An effective wing lift-curve slope , ap- 
plicable to basic lift distributions due to built-in twist, 
aileron deflection, roll, or aeroelastic twist, is approximately 
given at subsonic speeds by the relation 




A cos A 


A-p4-^ cos A 


( 1 ) 


where Ct^ is approximately given by 




1 — A// COS” A 


( 2 ) 


The basis of equation (1) is explained in reference 2. At 
supersonic speeds (more specifically, for supersonic leading 
and trailing edges), the effective wing lift-curve slope is 
approximately 


Cr._ - 


4 cos A 
■\ Mo^ COS" A — 1 


(3) 


provided A/<, is greater than 1/cos A. If A/^ is greater than 
1 but less than and not too close to 1/cos A, equations (1) and 
(2) may be used in the absence of better information; however, 
the results obtained for this range of Mach numbers should 
be used with caution. 

The lift-curve slopes given by equations (1) and (3) should 
not be confused with the rigid-wing lift-curve slope or the 
damping-in-roll derivative; they are merely effective values 
suitable for aeroelastic calculations. Values of the rigid- 
wing lift-curve slope and the damping-in-roll derivative can 
be used in conjunction with the methods of the present 
report and of reference 2 to obtain the values of the flexible- 
wing lift-curve slope and the damping-in-roll derivative 
because means are presented herein and in reference 2 for 
estimating the ratios of the flexible-wing to the rigid-wing 
values. For this purpose any experimental information 
concerning the rigid-wing values can be used; if none is 
available, references 3, 4, and 5 may bo used at supersonic 
speeds and reference 6, at subsonic speeds. 

The local aerodynamic centers are assumed to be at a 
constant fraction of the chord from the leading edge, so that 
thcA^ are all equal to the wing aerod 3 mainic center as a fraction 


of the mean aerod^mamic chord. The moment arm is 

then given b}” the relation 

= a (4) 

The local centers of pressure of the lift due to aileron de- 
flection are also assumed to be a constant fraction of the 
chord from the leading edge; and the moment arm €2 is then 
given by 

e2=cph—e (5) 

Theoretical two-dimensional values of the parameter cps arc 

presented in figure 2 for both trailing-edge and leading-edge 
ailerons at subsonic and supersonic speeds. At subsonic 
speeds the effect of finite span is to shift the center of pressure 
rearward. An appropriate value for this rearward shift 
may be estimated from the following relation, based on 
lifting-line theoiy for unswept elliptic wings with full-span 
ailerons: 


where the subscripts II and III refer, respect ivety, to two- 
and three-dimensional values. The use of the swept-span 
aspect ratio A\ in place of A should serve to extend this ap- 
proximate relation to swept wings. 

Theoretical two-dimensional values of aa, tlie angle of 
attack equivalent to unit aileron deflection, are also given 
in figure 2. At low aspect ratios the values of as for sub- 
sonic speeds tend to be higher than these two-dimensional 
values; as the aspect ratio approaches 0, as approaches 1, at 
least in the case of wings without reentrant trailing edges. 
Experimental values of both as and cps are preferable to 
theoretical values if the^' are available. For spoilers, ex- 
perimentallj' obtained values have to be used. 

The eftective lift-curve slope and the values of the section 
moment arms vary with the free-streain Mach number; 
hence, the appropriate values must be used at each flight 
condition for which aeroelastic calculations are made. 

The airspeeds at which the lateral-control aeroelastic 
phenomena are of interest enter the calculations in the form 
of the corresponding d 3 ’namic pressures. These d 3 Miamic 
pressures, in turn, are expressed in dimensionless form by 
means of the relations 



cos A 

144 

{GJ), 


Gl„ CiCrS,- COS A 

144 

{GJ)r 

<1 

Cl„ sin A 

144 

{El)r 


( 6 ) 

(7) 

( 8 ) 


The ratios of these quantities. 


^-q*~e,Cr{EI)A^'^ 

(9) 

fJ): tan A 
e eg C2Cr {EI)r 

(10) 


are independent of the d 3 Uiamic pressure and are veiy useful 
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for fumlyzing llio aoroelastic* boliavior of swept wings. 
olluM* dimoiisionloss panimelers, which are indepeudent of 
(he dynamic pressure, enter the calculations: 


and 



tail- A 

{ji'l Jr 


( 11 ) 

( 12 ) 


Structural parameters, — For the purposes of an aeroelastic 
analysis the wing structure is characterized b}' the location 
of the elastic axis, the magnitude and distribution of the 
bending and torsional stiffnesses El and (7«/, and the magni- 
tude of tlie rigid-body rotations imparted to the wing by its 
root (taken into account in this report only b}^ the location 
of an effective root). Tlie selection of these structural 
parameters is discussed in reference 2. 

PREUMINARV SURVEY OF LOSS IN LATERAL CONTROL 

'Pile information contained in some of the charts and ap- 
proximate formulas presented in the following sections of 
this report has been summarized in figures 3 to 5 for the 
purpose of ascertaining in advance of more detailed esti- 
mates, if desired, whether the aeroelastic phenomena con- 
sidered in this report are likely to affect the design of the 
wdng structure. This jireliminaiy suvve}' is not essential 
to an}" of the further calculations but may show them to be 
unnecessary in some cases. These figures pertain to constant- 
stress wings with half-span outboard ailerons. 

The charts of figures 3 (a), 4 (a), and 5 (a) pertain to wings 
of taper ratios 0.2, 0.5, and 1.0, respectively, with the moment 
arm co equal to 0 (corresponding roughly to subsonic flow 


conditions and elastic-axis locations faii*ly far back on the 
wing). These figures show" the dynamic-pressure parameter 
defined by either equation (6) or 




(1+X)2 




cos A 


18432 G n W hr 


Fs S 


Cr 


Fr1)a 


(13) 


plotted against tlie s^veep parameter k defined by erpiatioii 
(9) or 


l+XG Ax. . 

0 rp f tan A 

1 E Tji^ei 


(14) 


for several values of the aileron effectiveness parameter 
Ci^jCi^^ and for two values of the stiffness parameter 


or 


d_{EI)r /eiCrV 
\GJ)r \ S, ) 

d _ 32 ^ 

k- (1+X)- <? V 


( 15 ) 

( 16 ) 


Avhcl’e Fr is a root-stiffness parameter, and ? 7 a, rjbj Vsy Vsy 
Vq} Vioy n,nd t;i 9 are structural-effectiveness factors defined in 
reference 2. (If, at the time a preliminary survey of aero- 
elastic effects is to be made, no information whatever con- 
cerning the Aving stiffness is available, eqs. (13), (14), and 
(16) may be used; if an estimate of the root stiffnesses 
{GJ)r and {EI)r is available, eqs. (6), (9), and (15) may be 
used.) 

The charts of figures 3 (b), 4 (b), and 5(b) are the same as 
those of figures 3 (a), 4 (a), and 5 (a), respectively, except 
that these charts are for wrings with the moment-arm ratio 
6 equal to unity (corresponding rouglily to subsonic conditions 
wdth the elastic axis fairly far forward on the W"iiig). 



(a) Trailing-edge ailerons. (b) Leading-edge ailerons. 

Figure 2. — Theoretical two-dimensional values of the aileron-force parameters. 
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(a) Wings witli moment arm e-> = 0. 

Fkicrk o.— Charts for a preliminary survey of lateral control for flexible wings of taper ratio 1.0. 



Sweep parameter, k 

(b) Wings with moment-arm ratio €=1.0. 
Figure o. — Continued. 
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Ficiltre o. — Concluded. 


The charts of figures 3 (c), 4 (c), and 5 (c) pertain to wings 


with the moment arm Ci equal to zero (corresponding roughl 3 ' 
to supersonic flow conditions). These figures show the 
d 3 mamic-pressure parameter eg* defined 1 ) 3 ' either equation 
(7) or 

(^ ' Cl Co cos a 

. a) ^ 

^ 18432 G nW K j. 

Fu S Cr 

(17) 

plotted against the sweep parameter kje defined bv equation 
(10) or 

7 = - 2 - ^4. -—tan A 

(18) 

for two values of the stiffness parameter 


d {EI)r /e,tvV 
{kjef (GJ)r \ S, ) 

(19) 

or 


32 E 7767/7“ 

(20) 

{kje)- (1 +X)“ G \ Ax ) 77s77o77i5T7^>‘ 


The various lines of the charts of figures 3 to 5 designate 
the conditions at which a whig designed on ihe basis oj 
strength considerations alone is likely to encounter changes in 
aileron I’olling moment by various amounts due to wing 
flexibility. The lines for zero j-olling moment designate 
the conditions at aileron revei’sal. These charts should be 
used in conjunction with the preliminary survey charts in 
reference 2. The significance of tlie four quadrants of 
figures 3 (a), 3 (b), 4 (a), 4 (b), 5 (n), and 5 (b) is the same 
as that of the four quadrants of figures 2 (a), 2 (b), 2 (c), 
and 2 (d) of reference 2 and is discussed in the section “Pre- 
liminary Survey of Aeroelastic Behavior” of reference 2. 
The significance of the quadrants of figures 3 (c), 4 (c), 
and 5 (c) of the present report can be analyzed in the same 
wa 3 ^, except that the moment arm e-y takes the place of the 
moment Ci. These three figures are analogous to figure 2 (o) 
of reference 2 in that they pertain to the case Ci = 0. (Al- 
though k is infinite when c, is zero, which happens when the 
section aerod 3 ’namic centers are on the elastic axis, the 
parameter kje is not infinite in that case, except if €2 also 
happens to be zero, a condition which can be realized onl>" 
with a full-chord aileron.) 
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Aft or it lias boon foiiiul, through the use of the charts in 
reference 2, to what extent the wing design is affected b}’ 
such aeroelastic plienomena as divergence and aerod}'nainic- 
center shift, the same procedure can be used with the pre- 
liminary survey charts in this report to ascertain whether 
the wing design is likely to be affected by lateral-control 
difficulties. If these charts indicate the likelihood of sig- 
nificant aeroelastic effects on the aileron rolling moment, 
further calculations are desirable. The charts and approxi- 
mate formulas of this report ma}’ be used for the preliminary 
calculations; once the structure has been designed, more 
refined methods, such us that of reference 1, may be used. 

<'ALCULATION OF THE AEROELASTIC PHENOMENA RELATED TO LATERAL 

CONTROL 

Analysis of the many solutions for the aeroelastic phe- 
nomena considered in this report obtained by the methods 
given in the appendix shows that the data can be sum- 
marized by means of approximate formulas. These for- 
mulas involve the aerodynamic, geometric, and structural 
|)arametors of the wing through the dimensionless param- 
eters kj kje^ e, and d (eqs. (9), (10), (11), and (12), respec- 
tively) and through a series of constants Ki to K-. The 
constants are functions of the taper ratio, stiffness distribu- 
tion, and aileron span and are given in table I. As in 
reference 2, the form of these approximate formulas has 
b(‘en guided by considerations based on an idealized semi- 
rigid wing, and the actual values of the constants to 
K-; {Kt and K-y have been given in ref. 2) were obtained b}'' 
fitting the solution for the functions Bi to Bq defined in the 
appendix b}^ equations (AIS) and (A27) to their approximate 
(‘xpressions, equations (A32). 

Dynamic pressure at aileron reversal. — The solutions for 
the aileron reversal speed obtained by the methods given in 
the appendix can be summarized by approximate formulas 
which give the dimensionless parameters q*Rj («(Z*)«j or 
7/ji — that is, the values of the parameters defined in equa- 
tions (6), (7), and (8) which correspond to the value of the 
dynamic pressure at aileron reversal — in terms of the 
parameters kj kje^ 6, and d defined b}' equations (9), (10), 
(11), and (12), respectively. 

An approximate formula for is 


Tr~- 


K,(l-K,id) 


1 + {K, + K-Jud) 6 + K,d + K, y/ + K,h 


( 21 ) 


For veiy small values of the section moment arm (i and 
the resulting large values of the parameters e and k\ the 
following alternative forms of equation (21) are more 
convenient to use: 


^ 1 1 1 6 k 

-+K, + K2K,d^Ks - d + K, - j - 

€ e € K e 

and 

g,(i-g.f ./) 

\+{IU+K,K,d) l+Z-Ce 1 


( 22 ) 


(23) 


When the angle of sweep is zero, equations (21) and (22) 
reduce, respectively, to 


and 


* Kx 

(24) 

(«9 )h~ 1 

A%+- 

(25) 


and when the moment arm e\ is zero, as it ma}^ be in super- 
sonic flow, equations (22) and (23) reduce to 


and 


kA\-kAji) 

{eq ) r — — 7 

lU+K.Kid^K, - 

e 


qn— 




(26) 


(27) 


The constants K\ to Ki are given in table I for wings 
having half-span outboard ailerons and taper ratios of 0.2. 
0.5, and 1.0, for the two dift'erent types of stiffness distribu- 
tions. 33iese constants were found from the results of the 
numerical matrix method derived in the appendix. Also 
given in (able I are tliese constants for uniform wings liaving 
full-span ailerons and tip ailerons calculated from the results 
of the analytic integration method of the appendix. Since 
values of the constants K\ to AT; are given for three aileron 


TABLE I.— VALUJ-IS OF THE COEFFICIENTS Kx TO Ki 
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1 
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1 
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2.5S 1 

.381 
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' .972 

.715 1 

.194 

.028 

1.0 
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1.0 1 


1 1.0 
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.5 1 
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..T 

1.0 

1.2S7 
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' .410 
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. 2 
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! 
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spans in tlie case of the uniform wing, tliey may be inter- 
polated to 3”ield values for other aileron spans. No calcu- 
lations were made for other tlian lialf-span ailerons on 
tapered wings; nevertheless, as pointed out previously, qii 
calculated for half-span ailerons should be reasonably valid 
for outboard ailerons having spans wiiich differ considcrabl}'" 
from one-half. 

No aileron-reversal calculations have been made for swept 
wings with inboard ailerons. However, for unswept uni- 
form wings the dynamic pressure at aileron reversal lias been 
calculated for the limiting case of a wing with an inboard 
aileron of vanishingly small span by operational methods 
similar to those described in the appendix m connection 
with the calculations for uniform wings. The value of q*R 
obtained in this mamier is shown as a function of e in figure 
6. Also shown in figure 6 are the values of q*R for full-span 
and for tip ailerons. For small values of e, such as are likely 
to be encountered in subsonic flight, there is little difference 
m the values of q*R for the tlu*ee aileron configurations, but 
at large values of e, such as are likely to be encountered at 
high-supersonic speeds, there is some difference between 
them; the aileron reversal speed is highest for the tip aileron 
and lowest for the inboard aileron. However, these con- 
clusions ma}’ not be valid for nonuniform or swept wings. 

With the values of q*R, {eq*)R, and given b}' equations 
(21), (22), and (23) and the definitions of these parameters 
given by equations (6), (7), and (S), the values of q at 
aileron reversal may be determined. If desired, the corre- 
sponding airspeed ma}' be determined from the relation 

(28) 

If both qo (as obtained from ref. 2, for instance) and qR 
are positive and qR is greater than q^y there is no actual 
aileron reversal speed; in fact, tlie rolling moment due to 
aileron deflection will tend to increase with dAniamic pressure 
until divergence is reached. 




^ I 


. -2 




-4 



r 












IT 

1 










An.. 

Iwfl 






Divergence 









— 

Eq.(A 

36 ) 







-Full- 

Span aileron (s 
aileron (5/^—1 

■'= 0 ) 















-Tif 










J 






_ 















Inboard aileron 

^/"O, V-0)- 













■ — . 











n 





































k\ 













c. 

.11 





/ c 

"=C 

rN 














t:.. .-:i / . l\ 

} 

u- 

y 

\ 












111 

jnt;i 

Ull 






S 

jT 

'In 

_| 





(c. 

= 0 

LI 

, 5; 

[1 

« 

LJ 

• 0 ) 

— 









r 


11 

( 








! 

L_ 

j 


□ 

_J Li 

□ 

□ 

□ 

! 



-12 -10 -8 -6 -4 -2 0 2 4 6 

Moment -orm ratio, € 


8 10 


Fkiurk G. — Coinj)ari.<on of dynamic pressures at aileron reversal 
calculated by the analytic integration method of the appendi.x 
with those calculated from equation (A36) for unswept uniform 
wings. 


Nor is tliere an actual reversal speed if qR is negative, 
regardless of whether q^ is positive or negative, although the 
rolling moment due to aileron deflection will decrease slightly 
witli increasing d^^^amic pressure if qRjqo is greater than 1. 
Inasmuch as aileron reversal, unlike divergence, is not an 
instal)ilitv prolilem, although it is often convenient to anal^^ze 
it as such, there is no possibility of encountering aileron 
reversal in a liigher mode wlien qR is negative, at least not 
in tlie case of ordinary wings with straight leading and 
trailing edges and witli substantially straight elastic axes. 

The value of qn calculated for any given value of q*R, 
(€<2*)/e, or ~qR depends on the value of the effective lift-curve 
slope L\ and, hence, on the Mach number. As suggested 

in reference 1, the value of qR may lie plotted against Mach 
number on log-log coordinates; if the straight lines of the 
actual d^uiamic pressure at several altitudes as functions of 
Mach number are drawn on the same plot, an intersection 
of the reversal line with one of the lines of actual dynamic 
pressure designates possible aileron reversal at that value of 
dynamic pressure, Mach number, and altitude, 

Spanwise angle-of-attack distribution.— In the appendix, 
an approximate expression is determined for the change in 
angle of attack due to the deflection of ailerons on flexible 
wings. The ratio of the angle-of-attack distribituon due 
to structural deformation to the effective angle of attack 
of an aileron is 


^^qjq^ (K4 e + Kik) fa — Kjk Aj a 

j q 1 — Kok 

(Id 


(29) 


The functions /a and of the spanwise coordinate s* are 
given in figure 7 for wmgs having half-semispan outboard 
ailerons, taper ratios of 0.2, 0.5, and 1.0, and the two different 
types of stiffness distribution. The value of qo required in 
equation (29) ma}' be found from the approximate expression 
for q*o or q^ given in reference 2 as 


or 


2%= 


-K^k 


lu 



(30) 

(31) 


When qD is very large, a more convenient form of equation 
(29) is as follows: 


as ^q V 6 / 6 

(Xbb j ^ K\ 

qn 


(32) 


Spanwise lift distribution. — Within the limitations of tlic 
modified strip theory used in the analysis, the lift per inch 
of span is proportional to the local effective angle of attack, 
so that 


as 5 


=^qc 



(33) 
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whore aja^d is obtained as indicated in the preceding section 
and la is a unit step function of the distance along the span 
defined by 

2a=0 (when 5<s,) 

la= 1 (when 

Inasmuch as the functions /„ and A/^ have been calculated 
only for lialf-semispan outboard ailerons or spoilers, the 
expressions for angle-of-attack distributions and lift distri- 
butions (eqs. (29) and (33)) can be used diiTctly only for 
this case. 

Rolling-moment coefficient and rate of roll due to aileron 
deflection. — The rolling-moment coefficient due to a unit 
aileron deflection may be obtained in terms of its equiva- 
lent rigid-wing value from the approximate formula 




1-^ 
u 


C. 


1 --^ 

qo 


(34) 


and the wing-tip helix angle due to an aileron deflection 
pbj2Vy which is a measure of the rate of roll and the rolling 
maneuverability, may be obtained from 


7>6 

W 




Fkjl'ric 7. — Tlie angle-of-attack distribution functions and A/a for 
aileron deflections. s,* = 0.o. 



(b) Stiffnesses related to those given by constant-stress criterion for 
wings of taper ratio 0.5. 

Figure 7. — Continued. 



Dimensionless distance olong spon, s* 
(c) Stiffnesses proportional to cK 


Figure 7. — Concluded. 
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The manner in which the ailoron-efl'ectivoness parameter 
CiJCi^ varies with chTiamic pressure depends on the ratio 
of the reversal to tlie divergence dynamic pressure as may 
be seen from equation (34) and figure 8 . Wien the aero- 
dynamic center is ahead of the elastic axis, as is generally 
the case at subsonic speeds, qRlqo is positive and greater 
than one for swept forward wings, positive and less than one 
for unswejit wings, and negative for sweptback wings. 
Thus, in general, increases with dynamic pressure 

until divergence is reached for sweptforward ^\'ings; it de- 
creases slowly at first and then more rapidly as reversal is 
approached for unswept wings; and it decreases rapidly at 
first and then more slowly for sweptback wings. The rapid 
decrease in rolling eftectivcness for swept liack wings can be 
alleviated by the use of unconventional lateral-control 
devices which liave their centers of pressure ahead of the 
elastic axis, such as leading-edge ailerons or spoilers. These 
devices may also serve to make the dynamic pressure at 
aileron reversal negative, so that there is no reversal of 
lateral control in the given speed range. The loss or gaifi 
of lateral control is then given by the part of figure S for 
jiegative values of qhjR. 

The analysis summarized by the appi*oximate formulas 
(34) and (35) is based on rolling moments about the wing 
root instead of about the fuselage center line, partly to 
simplify the analysis and partly to avoid the introduction of 
the fuselage width as another independent parameter. How- 
ever, these approximate formulas should be valid for ob- 
taining rolling moments and rates of roll aliout the fuselage 
center line as well, because equations (34) and (35) are ex- 
pressed as ratios of flexible-wing to rigid-wing values; that 
is, the ratio of the rolling moment about the fuselage center 
line to its rigid-wing value should be nearly the same as the 
ratio of the equivalent rolling moment about the wing root 
to its rigid-wing value when the ratio of tlie fuselage width 
to the wing span is small. 

The rolling-moment coefficient and rates of roll given l)y 
equations (34) and (35) are functions of the aileron span 
inasmuch as qR is a function of the aileron span. The varia- 
tion of Ci^ or 2 ? 6 / 2 Hwith aileron span can therefore bo found 
onl 3 ^ for uniform wings; however, since the effect of aileron 



Dynamic -pressure ratio, <7/?^ 

FinuRB 8. — Variation of aileron rolling; effectiveness with ch'namic> 
pressure ratio and the parameter qnlqo- 


span on qR is not verv great, its effect on (\^ or pbl2V is not 
likely' to be great. The rolling moment due to the deflection 
of an inboard aileron can be found by superposition, because 
the aeroelastic ecpiations upon whicli the results of this re- 
port are based are linear (that is, the rolling-moment coeffi- 
cient due to a 30-percent -span inboard aileron, for example, 
is equal to (\^ for a full-span aileron minus Ci^ for a 70- 
percent -sjian outboard aileron), but onh’ in the case of the 
uniform wing is the required information for tlie full-span 
aileron presented herein. That the aeroelastic effecls on 
lateral control are likel\^ to be similar for inboard as for out- 
board ailerons, at least for unswept wings, ma}" be deduced 
from figure 6 . 

Inertia effects. — In steadv rolling flight no inertia effects 
are present which can affect the static aeroelastic problem 
except, possibly, for centrifugal forces on lieavy undcrslung 
nacelles, flflie maximum value of p 6 / 2 t^is therefore usuallv 
unaffected by inertia effects. However, in the equalh' im- 
portant problem of initial rolling acceleration, which governs 
the time in which a given rolling velocity can be attained, 
inertia effects must usually be taken into account. In 
reference 2 the observation was made that in symmetric 
flight inertia effects are not as important as other static 
aeroelastic effects, except for flving wings, because the inertia 
forces are in about the same ratio to the aero(h*namic forces 
as the wing weight is to the airplane weight. B}’ the same 
reasoning, however, inertia effects are almost alwaj’s verv 
important in getting into a roll, because the inertia forces 
are then in about the same ratio to the aerod 3 mamic forces 
as the moment of inertia of the wings about the longitudinal 
axis of the airplane is to the moment of inertia of the entire 
airplane about its longitudinal axis, a ratio which is usualH 
not much less than 1 . 

No charts are jiresented in this report for these inertia 
effects because the manner in which mass is distributed 
varies so wideh" among different airplanes that preparation 
of a generallA' applicable set of charts appears to be im- 
practical at present. However, the procedure outlined in 
reference 2 for taking inertia efiects into account in the cal- 
culation of quasi-static aeroelastic phenomena b 3 ' means of 
the charts presented therein ma 3 ’ be applied to the calcula- 
tion of the inertia effects encountered in starting a roll. 
This procedure is described in the following paragraphs. 

For a given rolling acceleration p, the linear normal accel- 
eration of an element of mass at a distance y from the center 
line of the airplane is py. From this linear acceleration and 
tlie known or estimated mass distribution of the wing, the 
inertia load /,- per incli of span and the inertia torque U per 
inch of span can be calculated for ain’ given normal, pitch- 
ing, or rolling acceleration. Substitution of these loads and 
torques for the terms / and le^c in equations (A3) or (A36) 
and equations (A 2 ) or (A35) of reference 2 , respect ivelA^, 
yields the values of the accumulated bending moments and 
torques in equations (A4) and (A5) or in equations (A37) 
and (A38) of reference 2 . Equation (A 6 ) of reference 2 , or 
the matrix equivalent of this equation, then 3 delds the anglc- 
of-attack distribution due to the deformations caused b 3 ^ the 
inertia effects associated with the given acceleration. 
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'Phis auglo-of-Httadv dislributiou can bo considered as a 
«:eo metrical aiigle-of-attack distribution; for the purpose of 
calculating tlie increinenl caused by aeroelastic action, this 
distribution can be approximated by a linear-twist angle-of- 
atlack distribution with a value at the wing tip which is such 
tliat the moment about the effective wing root of the area 
under the linear-twist distribution equals the moment of the 
area under the calculated aiigle-of-attack distribution due to 
inertia effects. (Tlie moment, rather than the area, is 
suggested as a basis of correlation because tlie angles of 
attack near the wing lip are more important in aeroelastic 
plienomena than those at tlxe wing root.) The justification 
for this rather arbitrary aixproximation to the angle-of- 
attack disti-ibution is that the correction to be applied as a 
residt of aeroelastic action to the deformations due to inertia 
loads is usually small compared Avith these deformations. 

1 'he angle of attack due to structural deformation 0:5 
associated with the linear-twist distribution can then be 
obtained from equation (21) and figure 7 of reference 2 or, 
if X — 0, from figure 8 of reference 2. The lift distribution 
associated with the total angle-of-attack distribution due 
to tlie deformations caused by the inertia effects, including 
tlie increment in this angle-of-attack distribution produced 
by aeroelastic action, (*an then be found from equation 
(24b) of reference 2 , in which ag and /o pertain to the calcu- 
lat(‘d angle-of-attack distribution due to the inertia effects 
(not the linear approximation to this distribution). This 
lift distribution can be integrated to obtain the rolling 
moment due to inertia effects, as modified hy aeroelastic 
action. 

fi'he rolling moment calculated in this manner may then 
be combined with the rolling moment due to aileron or 
spoiler deflection, which may be calculated as indicated in 
tlie ])receding section. If the contributions of the tail and 
tlie fuselage to the rolling moment are neglected. 


;)(/ — 7 ^)=^ Ci^^dqSb-\- 



wher(‘ f) is the angular acceleration in roll, / is the mass 
moment of inertia of the entire airplane about its longitudinal 
axis, is the mass moment of inertia of both wings about 
the longitudinal axis of the airplane, and (?/. is the flexible- 
wing value of the rolling-moment coefficient duo to aileron 
deflection (which may be calculated in the manner de- 
scribed in the preceding section). The ratio (d//'ir/c)p)s is 
the rolling moment per unit rolling acceleration due to 
inertia eflects, including aeroelastic eftects, and is equal to 
— plus the rolling moment due to the lift which results 
from tlie deformations due to the inertia loads per unit 
angular acceleration in roll as well as from the aeroelastic 
deformations wliich accompany these inertia deforma- 
tions; in other words, {dTJjcl^p)s is equal to —I,c phis the 
rolling moment calculated as described in the preceding 
paragraphs for /)=!. Then 
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(36a) 


or 


where 
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(36b) 


is a rolling-moment coefficient per unit aileron deflection 
which includes static aeroelastic effects, inertia eftects, and 
the aeroelastic magnification of the inertia effects. This 
rolling-moment coefficient is a truer index of the rate at 
which a roll can be initiated than Ci. . The initial rolling 

Og - 

acceleration (disregarding unsteady-lift effects) can be 
calculated from equations (36). 

illustrative example 

The approximate formulas described in the preceding 
sections have been used to find the effects of aeroelasticity 
on some lateral-control properties of the wing considered in 
the illustrative example of reference 2 , The resulting cal- 
culations are an extension of those in reference 2 , and the 
additional parameters are presented in table II. 


TABLK II,— PARAMETERS OF EXAMPLE \VIX(; 



1 Parameter 

1 

Subsonic (.\/= 0 ) 

Supersonic 

(.\/=1.5) 

■Si * 

0- 5 

0.5 

1 Si}* 

1.0 

I.O 

\ e>. . 

O.OIC 

0. 458 

- 

0. 0S2 

24. 1 

j k. .. 

7. 70 

79-0 


0. 551 

0. 551 

! A', 

2.82 

2.82 

: A-j 

0. 474 

0.474 

[ A '3 

0.917 

0.917 

A% 

1.000 

1.006 

A’s.. 

0.650 

0.650 

I\6. - 

0. IGO 

0. 160 

1 A-: 

0. 020 

0 . 0*20 

1 

-1.053 

-0. 0774 

^ K - 

LOST 

0.0713 

Ib/sq ft 1 

10. 300 

2.500 


Sec fiR. 7(c) 

See fig. 7(c) 

AA 

See fig. 7(c) 

See fig. 7(c) 

1 a, laid 

Ct,/c,, - - 

* *0 

See cq. (29) 
l-f-0.024 

QP 

i-i 

Qd 

.See c(i. (29) 

1+1.085 — 

1 -^ 

I 90 

; (pbl2V) 

1+0 024 — 

1 + 1.0S5 — 

cpW 2 r)o 

7/> 

flit 
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The subsonic and supersonic values of the parameters 
6, and d were calculated from equations (9), (11), and (12), 
respectively. With the appropriate values of the factors Kx 
to K-j interpolated fx-om table 1, the values of <i^r wei*e cal- 
culated from equation (21) and are given in table II. From 
these values of the subsonic and supersonic dynamic 
pressures at aileron revex'sal wei*e found b}’ means of equation 
(6). These values of ^r vaiy as the reciprocal of the effective 
lift-curve slope, if the corresponding values of and are 
assumed to remam constant. 

In order to find the angle-of-attack distributions due to 
deflections of the aileron from equation (29), the values of 
the functions /a and A/„ were taken from figure 7(c). The 
spanwise change in angle of attack is shown in the top plot 
of figure 9 for different values of the d^mamic-pressure ratio. 
The rolling-moment coefficient and wing-tip helix angle due 
to deflections of the aileron were calculated from equations 
(34) and (35) and were plotted in figure 9 as functions of the 
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Figi'Re 9. — Effect of aeroDlastic action on some lateral-control prop- 
ortie.s of the example wing. 


DISCUSSION 

LIMITATIONS OF THE CHARTS AND APPROXIMATE FORMULAS 

The charts and the approximate formulas presented in 
this report arc subject to certain limitations as a result of 
the approximations made in the calculations on which they 
are based. Those limitations are discussed fully in reference 
2 and can bo classified as restrictions on the plan form, on 
the speed regime, and on the wing struct lire. The limitations 
are given very briefly as follows: 

(1) The results obtainable by the use of the charts and 
approximate formulas are likely to be unsatisfactory for 
wings of very low aspect ratio, veiy large sweep, or zero 
taper ratio. 

(2) The results are restricted to wings on which tlie span- 
wise lift distribution is roughly proportional to the chord 
and angle of attack and on which the section aerodynamic 
centers are at an approximately constant fraction of the 
chord; these restrictions are most likely to be violated by 
wings flying at transonic speeds and by wings liaviug con- 
centrated sources of lift, such as nacelles and tip tanks. 

(3) Tlie results are somewliat resti-icted to wings with one 
of the two types of spanwise stiffness distributions used 
in the analysis, wings with no chordwise bending (relatively 
thick wings), and wings having an elastic axis at an approxi- 
mately constant fraction of the chord. 

The manner in wliich the aeroelastic effects of aileron 
deflection arc analyzed in this report imposes certain ad- 
ditional limitations that particularly affect tlie aileron 
geometry. In the analysis in the appendix, the spanwise 
lift distribution due to the deflectioii of an aileron was 
approximated by strip theory, an approximation which is 
probabh’ less valid for aileron deffectious tliaii for geometric 
angles of attack. The assumption was also made that the 
centers of pressure due to aileron deflections are at a constant 
fraction of the wing chord; this assumption is also probably 
less valid than the assumption that the section aerod\mamic 
centers of the wing are at a constant fraction of the chord. 
Since these assumptions are more nearly true for Avings of 
high than those of low aspect ratio, these limitations serve, 
m effect, to restrict the applicability of the px-esent report 
to aspect ratios somewhat higher tJian those amenaiile to 
the analyses of reference 2. 

The results of the present report do not take into account 
explicitly any flexibility of the aileron itself because of the 
assumption that the angle between the aileron and wing is 
constant along the span of the aileron. This assumption is 
almost universally made in anal 3 'zing the aeroelastic prop- 
erties of ailerons and is justifiable because the net effect of 
the difi’ereiice between the wing deformations and the aileron 
deformations on the overall lift and moments appears to 
be negligible. 

As a result of the fact that the static aeroelastic phenomena 
associated witli lateral control involve many more param- 
eters than do those associated with sA’mmetric flight, tlie 
coverage of the various parameters is not as complete as 
in reference 2. Spccificalh', all the charts and approximate 
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foriiuilas jirosoiiled herein arc restricted to outboard lateral- 
coiilrol devices (ailerons or spoilers) except for the uniform- 
wiiifi: case, and most of the calculations upon which these 
r(‘sulls arc based were made for wings with half-semispan 
outboard ailerons However, the results 

of the analysis for the uniform wings with full-span ailerons 
0s‘,* = 0, and tip ailerons indicate 

that, except for the angle-of-attack distributions, the results 
based on a half-span aileron are approximately valid for 
outboard ailemns of spans differing considerably from one- 
half. The ratio of the ailemn chord to the wing chord is 
assumed in the charts and formulas to be constant over the 
sj)an of the aileron. 

RELATION RETWEEN STRENGTH AND STIFFNESS AS DESIGN CRITERIA 

'I'he relation between strength and stiffness as design 
criteria was discussed in reference 2. The preliminary sur- 
vey (‘harts in reference 2 indicate the extent to which a 
wing designed on the basis of strength considerations alone 
is lik(dv to be affected by aeroelastic phenomena and, con- 
secjuently, indicate whether the wing has to be stiffened 
l)(*vond the amount associated with the required strength, 
ff'he jireliminary survey charts of this report (figs. 3 to 5) 
serve the same purpose in regai-d to the aeroelastic effects on 
lateral conti’ol; furthermore, even though the charts of 
reference 2 may indicate that a particular wing is not sig- 
nificantly affected by the aeroelastic phenomena considered 
in that report, this wing may still have to be stiffened be- 
cause of an undesirably large loss in lateral control. 

As may bo concluded from the survey charts of reference 2 
and I he present report, as well as from the discussion con- 
tained in reference 2, the following wings designed on the 
basis of stiTUgth considerations are most likely to be subject 
to adverse aeroelastic effects on lateral control and rolling 
maneuverability: 

(1) Wings operating at a high speed or dynamic pressure 

(2) Sweptback wings 

(3) ffdiin wings 

(4) Wings designed for low wing loading 

(5) Unswept and moderatel}^ swept wings with an elastic 
axis relatively far forward on the chord or with the center 
of pressure of the lift produced by aileron deflections rela- 
tively far back on the chord as a result of the aileron configu- 
ration (small aileron chord or wing of low aspect ratio) or 
flight condition (supersonic speeds) 

(6) Wings witli a relatively liigh lift-curve slope 

STRUCTURAL WEIGHT ASSOCIATED WITH THE REQUIRED STIFFNESS 

WHien a given wing has been shown to be subject to un- 
desirably large aeroelastic effects b.y means of the charts of 
reference 2 and this report or b}" any other method, the 
problem arises how to distribute the additional required 
stiffness — that is, which spanwise distribution of structural 
material will alleviate the adverse aeroelastic phenomena 
to the desired extent with the minimum increase in struc- 
tural weight. 

In order to shed some light on this problem, aeroelastic 


and weight calculations have been made for wings with a 
taper ratio of 0.5 with a family of somewhat arbitrarily 
selected stiffness distributions which differ from the distri- 
bution required by the constant-stress criterion in a manner 
described in reference 2. These stiffness distributions are 
designated by the tip stiffness ratio oj, which is the ratio of 
the stiffness El or GJ at the wing tip to the corresponding 
stiffness of a constant-stress wing. The results of the 
lateral-control calculations for wings with taper ratio 0.5, 
with constant wing-thickness ratio hic along the span, and 
with two of these stiffness distributions are included in 
table I and figure 7 (b). The designation ^^excess strength” 
refers to the stiffness distribution increased over the constant- 
stress recpiiremcnt to such an extent that the value of is 
2.0. The results of the aeroelastic calculations for the 
stiffness distributions decreased below the constant-stress 
requirement to a value of co=0.5 happen to be the same as 
the results for the constant-stress stiffness distributions for 
wings with linearly varA-ing wing-thickness ratio and 


The structural weight considered in these cal- 

{h!c)r 

culations is that of the primary load-carrying structure; 
the remaining structure is assumed to be unchanged in the 
stiffening process. 

The results of the weight calculations and aeroelastic 
calculations in reference 2 indicated that in tlie case of wings 
with taper ratio 0.5 stiffening the structure in the outboard 
region of the wing (co greater than 1) was more efficient, 
from weight considerations, in alleviating the aeroelastic 
effects considered in that report than the addition of stiffness 
in the inboard region. This conclusion is corroborated, in 
essence, b}’ the calculations made for the aeroelastic phe- 
nomena considered herein. Figure 10, which consists of a 
plot of the structural weights required for a given loss (20 
percent) in lateral control at a given dynamic pressure, indi- 
cates that the least weight is associated with values of the 
tip stiffness ratio co greater than 1, except for wings with 
values of the sweep parameters k or Ar/e equal to —8. These 
large negative values of k or kje^ however, pertain to wings 
that are either (1) sweptforward, or (2) sweptback with the 
aerodynamic center behind the elastic axis (in the case of 
negative Ar), or (3) sweptback with the center of pressure due 
to aileron deflection ahead of the elastic axis (in the case of 
negative A^/e, as it ma}^ be for spoilers or leading-edge 
ailerons) . 

For sweptforward wings, the aileron rolling moment 
usually increases rather than decreases with dynamic pres- 
sure, so that lateral control does not impose any structural 
requirements. For sweptback wings with negative values of 
Cl and Co, alleviation of the aeroelastic effects in lateral con- 
trol can be effected at the least cost in weight b}' adding 
structural material in the inboard region of the wing; inas- 
much as the minimization of the shift of the aerodynamic 
center of these wings can be effected most efficiently b}" 
stiffening the outer region of the wing, a compromise must 
be made if both t 3 ^pes of static aeroelastic phenomena are of 
concern. 



18 


REPORT 1139 NATIONAL ADVISORY COMMIl'TEE FOR AERONAUTICS 



(a) f>=0. 

(b) e=l. 

(0 f,=0. 


in reference 2 ^vas the fact that even if the conditions of 
equation (37) are achieved there may still be great losses in 
lateral control, and the wing ma}^ still be subject to adverse 
d^mamic phenomena; in fact, the severity of adverse aero- 
elastic effects in the lateral control and of certain dynamic 
phenomena may be increased as a result of achieviog aero- 
isoclinicism. 

The results of the present report corroborate the conclu- 
sion concerning aeroelastic effects on lateral control. How- 
ever, by suitable additional modifications a Aving which has 
been made aeroisoclinic can also bo made to suffer no loss in 
lateral control due to aeroelastic action. As may be seen 
from figures 3 to 5 or from equation (34), the condition for 
no losses in lateral control is that be equal to qo or, by 
setting equation (21) equal to cc^uation (30), that 

+ I ^/=0 (38) 

of which the condition 


{GJ)r , , 1 

‘ — ' tan A= 

€2Ct {EI)r K 2 


(39) 


usually is an approximate solution. Therefore, in order to 
satisfy both of the conditions specified in equations (37) and 
(38), the section moment arm €2 must be nearly equal to 
— Cl or, in other words, the center-of-pressure parameters a 
and cp 5 must be nearly equal, a condition Avhich can be 
satisfied b}' using full-chord aileroiis (all-movable wing tips) 
or a combination of geared leading- and trailing-edge ailerons, 
for instance. HoAvever, as pointed out in reference 2, 
attempts at soh’ing static aeroelastic problems by aiming at 
aeroisoclinicism may tend to aggravate certain dynamic 
phenomena. The same statement must also be made con- 
cerning the foregoing methods of alleviating static aeroelastic 
effects on lateral control; these methods may, for instance, 
lead to flutter difficulties Avhich may require excessive mass 
balancing. 


Figure 10.— The effect of tip stiffness ratio on tlie structural weight 
required to maintain a constant level of lateral-control effectiveness 
^Cfj = 0.8C4^) at a given dynamic pressure for wings with a taper 

ratio of 0.5. 


SOME REMARKS CONXERNING THE AEROISOCLINIC WING 

As shown in reference 2 an overall type of aeroisoclinicism, 
in Avhich bending and torsion action tend to cancel for the 
wing as a whole, can be achieved for the aeroelastic phe- 
nomena considered in that report by a choice of a suitable 
ratio of the bending to the torsion stiffness or by a choice of 
the elastic-axis location, that is, by satisfying the relation 


(GJ)r 
eiCr {EI)r 


tan 



(37) 


RELATION OF THE CHARTS TO DESIGN PROCEDURE 

The conventional procedure of designing a wing on the 
basis of strengtli requirements and later checking it for 
aeroelastic effects can be facilitated at several stages by 
using the methods described in the present report and in 
reference 2. As pointed out in reference 2, for instance, the 
preliminary -survey charts presented therein can be used to 
establish some static aeroelastic characteristics that Avould 
be obtained in symmetric flight if the Aving Avere designed for 
strength alone. 

If these characteristics are deemed satisfactory, the de- 
sign can proceed on the basis of strength requirements alone. 
If, on the other hand, they are considered unsatisfactory, 
the Aving must be stiffened. The amount of additional 
material required can be estimated, as indicated in reference 
2, by interpolating betAveen the results presented therein 


(AAdiere K 2 is given in table I). HoAvcA^er, also pointed out 
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for tlio constant-strength case, the “excess strength” case, 

and the case previously mentioned, the 

additional structural material is usuall}’' most effective if 
distributed near the wing tip. 

Similarlv, the preliminary-survey charts of this report 
can be used to ascertain whether the wing can meet lateral- 
control requirements if designed for strength alone. If it 
must be stiffened to meet these requirements, the necessaiy 
amount of additional material can be estimated in the 
same manner as indicated in reference 2 for aeroelastic 
(dfects iiK'urred in symmetric flight. 

Inasmuch as the chai*ts of reference 2 and of the present 
n‘port pertain only to static aeroelastic phenomena, the 
problem remains of ascertaining in the preliminary design 
stage whether a wing designed for strength alone (or, for 
(hat mattei*, a wing designed both on the basis of strength 
n^qiiij-ements and of static aeroelastic considerations) is 
likely to e.xperience flutter difficulties. However, flutter is 
a much more complicated phenomenon and depends on 
many moj*e parameters than do static aeroelastic phenomena. 
(\)nsequently, preparation of a generally applicable set of 
charts appears impractical. Nonetheless, although phe- 
nomcnologicall}^ flutter is not related to the static aero- 
(‘lastic phenomena considered in reference 2 and herein, it is 
mechanically related bv virtue of the fact tliat all these 
l)henomena depend on the wing geometiy and the wing 
slitfness (although the aerodynamic parameters are different 
and lluttei*, unlike the static aeroelastic phenomena, involves 
(he mass distribution of the wing and the damping properties 
of the structure). On the basis of past experience, certain 
qualitative cojiclusions can be drawn concerning this 
relation. 

As shown in the charts of this report, the aileron reversal 
speed, or the speed at wliich a specified amount of control is 
retained, is lower for highly sweptback wings than for 
imswept wings. Similarly, the divergence speed decreases 
rapidly as the angle of sweepforward increases. For a 
typical family of wings the values of at divergence and 
at reversal as obtained by the charts herein are shown as a 
function of the sweep parameter k in figure 11. For un- 
sw(‘pt wings the (Kmamic pressure at flutter is usually within 
a c(‘i*tain range varying between a value lower than the 
dynamic pressure at reversal to a value higher than the 
dynamic pressure at divergence, depending on the geo- 
metric, structural, aerodynamic, and mass parameters of 
(he given case, and varying even for a given case and a 
giv(^n speed range with altitude, because a change in air 
(l(‘nsity ma}" change the mode in which the wing flutters, 
fl'his wide range is iiulicated in figure 11 b}' starting three 
flutter curves (which do not necessarily describe the upper 
and lower limits) at k=0. If the variation of flutter speed 
with sweep angle is assumed for the purpose of illustration 
to b(* similar to that indicated in figure 17 of reference 7, 


the flutter curves of figure 11 are obtained. 

This figure must not be construed as presenting any 
quantitative information; to emphasize this point the family 
of wings is not identified. Even qualitative!}' the relation 
between the dpiamic pressures at flutter, divergence, and 
reversal is subject to certain limitations because the flutter 
tests of reference 7 were performed at subsonic speeds on 
models without ailerons and concentrated masses, which 
fluttered in the classical two-degi-ee-of-freedom mode. 1'hero 
is reason to believe that sweptback wings with high aspect 
ratios flying at high altitudes may experience a possibly 
mild form of flutter in a single-degi*ee-of-freedom mode, be- 
cause a vertical motion necessarily implies vertical bending 
and, hence, in the case of a swept whig, a variation in tlie 
angle of attack. In general, the gi-eater the number of 
degrees of freedom the more difficult it is to relate flutter to 
the static aeroelastic phenomena. 

However, at subsonic speeds and low or moderately high 
altitudes at least, the trend shovui in figure 1 1 sliould be valid 
for wings without very large concentrated masses and with 
inwersible controls, whicli tend to minimize the possibilit}’ 
of aileron-coupled flutter. Consequently, if these wings are 
liighly swept back they can be designed to meet lateral- 
control requirements with the likelihood that they will then 
be safe against flutter as well, provided conventional lateral- 
control devices are used. On the other hand, if these wings 
are unswept or even moderately swept back, they may have 
to bo stiffened beyond the amount required by static aero- 
elastic considerations, or mass balanced, or both. 

In any event the final desigti must be checked both for 
static aeroelastic effects and for flutter. In many cases the 
static aeroelastic effects can probably be calculated with 
sufficient accuracy by means of the charts and approximate 
formulas of reference 2 and the present report. In some 
cases, however, particularly if these effects are in any way 
critical, a more refined method of analysis, such as that of 
references 1 and S, may have to be used. 



Sweep porometer, k 

Figure 11. — Variation of critical dynamic pressures with the sweep 

parameter k. 
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CONCLUDING REMARKS 

An approximate method based on charts and approximate 
formulas has been presented for estimating rapidly the 
aeroelastic effects on the lateral control of swept and unswepfc 
wings at subsonic and supersonic speeds. The chaids and 
approximate formulas presented herein together with those 
presented in NACA Report 1140 also serve to simplify 
design procedure in many instances because the}^ can be 
used at the preliminary design stage to estimate the amount 
of additional material required to stiffen a wing which is 
strong enough and because the}' indicate that the best way 
of distributing this additional material in most cases is to 
locate most of it near the wing tip. 

For the purpose of making specific calculations, the limita- 
tions of the method of this report are that they do not apply 
dii'ectl}^ to wings with very low aspect ratio, witli very large 
angles of sweep, with zero taper ratio, or with largo sources 
of concentrated aerodmamic forces. 


The charts and approximate formulas indicate that the 
control effectiveness of an airplane may be increased l>y 
varying some of the design parameters such as the ratio of 
torsional to bending stiffness and, if necessary, resorting to 
unconventional lateral-control devices. The charts also in- 
dicate tliat a wing which is strong enough is most likely to 
be affected by losses in lateral control due to wing flexibility 
if it is to operate at high d}uiamic pressures, if it is thin, if 
it has a large angle of sweepback, if it has an elastic-axis 
location relatively far forward on the chord or a location of 
the center of pressure due to aileron deffection far rearward 
on the chord, or if it is to operate at transonic or high super- 
sonic Mach numbers. 


Lanoley Aeronautical Labok.vtorv, 

Xational Advisory Committee for Aeronautics, 
Lanoley Field, Va., Ma?‘ch 7, J952, 
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APPENDIX 

METHODS OF CALCULATIONS ON WHICH THE CHARTS ARE BASED 


THE AEROELASTIC EQUATIONS 

'I'hc assumptions made in the following analysis are the 
same as those made in reference 2: 

(1) Aerod^mamie induction is taken into account b}* 
applying an overall correction to strip theoiy. 

(2) Aerodjmamic and elastic forces are based upon the 
assumption of small deflections. 

(3) The wing is clamped at the root perpendicular to a 
straight elastic axis, and all deformations are considered to 
be given b}^ the elementaiy theories of bending and torsion 
about the elastic axis. In addition it is assumed, as in 
reference 1 and elsewhere, that the angle between the 
aileron and the wing is constant along the span of the 
aileron. 

flMien, for a wing with an outboard aileron, the force per 
unit width on sections perpendicular to the elastic axis is 

qcC\^ 

(aj + aaSio) (Al) 


where is a unit step function of s defined by 
la(s) 


H 


where Si is the s pan w ise ordina te of the inboard end o f the 
aileron. The running torque of this force about the elastic 
axis is 

— (as — € aid la) (A 2) 


where e is the moment-arm ratio eojei. 

The integration of these forces jdelds the accumulated 
tonpie and bending moment: 


qCr^e^Cr 


A'I= 


(ICtOl 


m 2 

(a. — < 


ai5la)ds 


as “h ai8l(^ (Is (Is 


(A3) 

(A4) 


Combining these equations with the equations of elastic 
deformation presented in appendix A of reference 2 as 


’~Lgj 


Tds 


M ds 


results in two simultaneous equations of equilibrium: 


A 

ds 



[(^ cos A— r sin A) — eagSifl] 


[(<p cos A — r sin A)-ro:j5ial 


(A5) 


(A6) 


These equations are subject to the following boundary 
conditions: 


^(0)=0 


(A7a) 

r(o)=o 


(A7b) 


= 0 

(A7c) 

'n 

II 

= 0 

(A7cI) 

II 

= 0 

1 

(A7e) 


The angle of attack due to structural deformations is 
related to <p and F by the equation 

as = (p cos A — r sin A (AS) 

After equations (A5) and (A6) have been solved, the rolling 
moment about the wing root may be found from the ex- 
pression 

qS Cl 

— ^ = Tr sin \+Mr cos A (A9) 


where the root twisting moment Tr and the root bending 
moment Air are given by equations (A3) and (A4) evaluated 
at s equal to zero. Then the rolling-moment ratio becomes 




^=1- 


i eBs+B, 


(AlO) 
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where the functions B\ to are defined b}’ 



(Alla) 

B,+eB,== f‘ r** y ds* ds* 

do d 0 0^50 

(Allb) 


(Allc) 

5a= r f '*- la ds* ds* 

dodo f' r 

(And) 

and where the parameter d/k is defined b}- 


d €\Cr . 
v== — tan A 

k St 

(A 12) 

SOLUTION FOR UNIFORM WINGS 


If the torsional stiffness, the bending stiffness, and the 
chord of the wing liave constant values of {GJ)r, {EI)rj and 
Cr, respectively, along the wing span, the equations of equilib- 
rium (A5) and (A6) become 

cos A = — g*[((^ cos A — r sin A) — €as5 id 

(A 13) 

sin A=“^ [(9? cos A— r sin A) + ai5 id 

(A14) 

where the differentiation denoted by the prime is 

with re- 

spect to J = 1 

S t 


Differentiating equation (A 13) once with respect to J and 
combining it with equation (A 14) 3’ields the single differential 
equation of equilibrium, 

?0's=a55 g*e i„'+<7 la 

(A15) 

subject to the following boundaiy conditions: 


«.(!) = 0 

(A 16a) 

a/(0) = 0 

(A 1Gb) 

a/'(0) ——q^ [«s(0) — e a^h la (0)] 

(A 16 c) 


and /3±/7. Those tlii’oe functions are equal to zero when 
KO- 

The substitution of this solution into equations (All) 
yields the funct'ons 

+ (AlSa) 
(AlSb) 

B.{cf, - f.)- ( AlSc) 

BM*, ^')=p [/3(1)-/3(1 -?,)1 + 

[p--{j[Jj][/5(l)-/3(l (Aisd) 

B~-^i (Al8e) 

B,= ^i-\ir (Also 

Tlie value of (/* at reversal is that value which makes 
in equation (A 10) equal to zero for ^iven values of 
the parameters k, €, and d. 

For a full-span aileron ($i— 1, Jo — 0), equations (A18) 
become 

(AlOa) 

/iA9^^0-J[.^(l)-ll+[|-■^;[jj]/3(l)<7* (AlOb) 

B,(q*,k)=-j^y^+B, (Al9c) 

BM*^ 1] + j(.A(l)- l]+[p-^J]]/3(l)T 

(AlOcl) 


The complete solution of equation (A 15) can be readily 
obtained by means of Laplace transforms as 

— ;^=|^/.-AJ)~/5(f— Jd — /s(J) J (l-r€)(7 -f 

/3©-/3(s^-fd-2„(|) (A17) 

where the functions f^{^) and /oCj), as well as /4(J) which 
will be used subsequently, are defined in aj)pendix A of 
reference 2 as 

/3(?)= ae-=«+e'*«(c', cos sin 7?) 

/i(?) = C^e -f ( ^5 cos 7 $+ y sin 7 ?) 

/5(?)=C:C--«-bc«(Cs cos 7?+^ sin 7?) 

where the constants of integration are defined in reference 
2 in terms of the roots of the characteristic equation, —2/3 


B,^-l (Al9c) 

^6=^ (Al9f) 

and for a “tip aileron” of very short span (Ji-^0, J<,= 0), 

equations (18) become 

/A(</^^0={/3(l) + [/4(l)-^7.(l)];^■-^] r/-l| ^,+B, (A20a) 

5,(g = [2 \f:,( 1) - 2* )] (A20b) 

+ (A20c) 

B,{q*,k)^ j '' 

B,= -^i (A20e) 

/?o= e. (A200 
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SOLUTION FOR NONUNIFORM WINGS 

By moans of strip thoorv applied at a finite uuniber of 
pointv, ocjiiation (Alj may bo written in matrix notation as 

(A21) 

and the expression for torque (A2) as 

(A22) 

1'lie matrix notation used in this analysis is the same as in 
reference I . 

E(|uations (A3) and (A4), written in matrix notation, arc 


{ = } (A23) 

{ } = T44 ' L^l { { H (A24) 

wh(‘i*e the integrating matrices [/'] and [//'] are, respec- 
tively, for single and double integration from tip to root 
and are given in reference 1. The value of {2^} at the 
matrix station nearest the discontinuity of la can be modi- 
fied as in reference 1 in such a manner that preinultiplication 
of \ la} by [7hJ and [Il\\ yields the same area and mo- 
im‘nt about the wing root, respectively, as would be obtained 
l)V analytic integration. This modification may a'so be re- 
garded as an attempt to round off the lift distribution near 
tlie inboaixl end of the aileron in a physically reasonable 
manner. 

ff'hc combination of equations (A23) and (A24) with 
matrix expressions for {<p} and {F} in terms of {T] and 
{ M] (see appendix A of ref. 2) yields the equilibrium equation 

[1 1 1 - ry*[A]] {«,} = - q*a,5 { ^ } (A25) 

which can Ijc solved for { aj}. The aeroelastic matrix [A] is 
defined in reference 1 and the column matrix {d} is defined 
hy 




-f '1 [^]] f ) »26) 

After the rolling-moment ratio expressed b}’ equation 
(A 10) is set equal to zero, the condition for aileron reversal 
is 

where the matrix equivalents of equations (All) become 



(A 27a) 


(A 2 7b) 


(A27c) 


B,^\II 


'■'Ul*''-’ 


(A27d> 


Therefore, the condition for aileron reversal expressed in 
matrix notation is 

Solving equation (A28) for and multiplying the resulting 
equation by {j8} yields 




The substitution of this expression for into the 

equilibrium equation (A25) yields 

{ 0 : 5 } = [ Aie] { Of, } ( A30> 

where (he aileron-reversal matrix [Ar] is defined by 
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1(0’ 
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(04 

-iii'A 

C~] 
pr _ 
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<g)[ri^'4(0'l+'"'44 

i‘‘''4(0T‘'''40J''-> 


(A31) 


The value of 2* tit aileron reversal can be found by the itera- 
tion of this aileron-revei'sal matrix. 

REPRESENTATION OF RESULTS BY APPROXIMATE FORMULAS 

Approximate formulas, similar to those in reference 2, 
have been used to combine the results of the many compu- 
tations indicated in this anah'sis; as in the case of the for- 
mulas presented in reference 2, those presented in this sec- 
tion are based on considerations of a semirigid wing. 

The functions Bi to in equation (AlO) have been found 
to be given quite accurately by the following approx mate 
expressions: 


B, = -B, 


-t * 




- ^^\_q*_^-K2k 

I * 


B,= -B, 


q ■ D 

q*lq*D {K.^IQk 
, q* l-K^k 


(A3 2a) 


(A 3 2b) 


(A3 2c) 


B,= -B, 


D 


q* l~IC 2 k 


q ' D 

B~-B,K, (A32e) 

wliere the factors to iiT: depend on the aileron span, the 
taper ratio, and the spanwise variation of the bending and 
torsion stiffnesses. The factor Kj is independent of the 
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aileron span, and the factor is independent of the stiffness 
variation. 

The following approximate formula may by obtained by 
the substitution of equations (A32) into equation (AlO): 


C, 


1-^ 
^ (jR 


1 


(Id 


where the value of t? at aileron reversal is 


(A33) 


The results of reference 2 indicate that the damping-ia-roll 
derivative may be expressed approximate!}' by 


Cl 


1 



qo 


and since the wing- tip helix angle due to roll is 


'ph 
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('k 


(A37) 


(A38) 


qn= 


qo(l-K,k)(^l-K,jd^ 


l+{K,+K,K,d)e+K,d+K, j d+K,k 


(A34) 


an approximate formula for rolling maneuverability is 


pb 

2F"' 




(A39) 


With the use of the approximate formula pr<‘sen ted in 
reference 2, 


* 

9 D 


K, 

i~K.dc 


(A35) 


equation (A34) becomes 


K 


q*R= 




1 +{K,+K J<::,(l)e~\~K,(l + K, ^ d+K,k 


(A36) 


The accuracy of equation (A36) compared with the results 
calculated directly by the method of the preceding section 
is illustrated in figures 6 and 12. 



Figure 12. — Comparison of dynamic pressures at aileron reversal 
calculated by the analytic integration method of the api)endix with 
those calculated from equation (A3G) for swept uniform wings with 
tip ailerons, s,*— 


Figure 13 shows the approximate formulas (A33) and (A39) 
to be in good agreement with more accurately computed 
values. (Actually, a slightly more accurate method of 
estimating Cj is given in ref. 2 but the simpler expression 
given here is comparable in accuracy to the expression for 



Figure 13. — Comparison of rolling power and rolling maneuverability 
calculated by the matrix method of the appendix with those cal- 
culated by equations (A33) and (A39) for uniform wings. s,* = 0.5; 

e=1.0. 
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Ci^, and its simplicity facilitates its application to the esti- 
mation of the rolling maneuverabilit 3 ^) 

An approximate expression for the structural twist due 
to aileron deflection similar to the expressions for the struc- 
tural twists due to geometrical angles of attack given in 
reference 2 has been deduced from the results of the anal}’sis 
in the preceding section: 

= g/gp -\-K-2k)ja — KokAfa ('A40) 

I ^ ^ — Kok 

Qd 

where /« and A/a are functions of the spanwise coordinate 
the wing chord and stiffness variations, and the aileron 
span. The accuracy of equation (A40) is indicated in 
figure 14. 
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Figure 14. — Comparison of angle-of-attack ratios calculated by the 
matrix method of the appendix with those calculated by equation 
(A40) for uniform wings. s,*=0.5; e=1.0. 
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